Electron relaxation in metals: Theory and exact analytical solutions 
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The non-equilibrium dynamics of electrons is of a great experimental and theoretical value provid- 
ing important microscopic parameters of the Coulomb and electron-phonon interactions in metals 
and other cold plasmas. Because of the mathematical complexity of collision integrals theories of 
electron relaxation often rely on the assumption that electrons are in a " quasi-equilibrium" (QE) 
with a time-dependent temperature, or on the numerical integration of the time-dependent Boltz- 
mann equation. We transform the integral Boltzmann equation to a partial differential Schrodinger- 
like equation with imaginary time in a one-dimensional "coordinate" space reciprocal to energy 
which allows for exact analytical solutions in both cases of electron-electron and electron-phonon 
relaxation. The exact relaxation rates are compared with the QE relaxation rates at high and low 
temperatures. 
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I. INTRODUCTION 

In recent years investigations of photo-response func- 
tions in advanced materials have gone through a vigorous 
revival. In particular, laser "pump-probe" techniques, 
where a second probe pulse is delayed in time with re- 
spect to the pump pulse, provide unique information 
on the strength of electron-electron (e-e) and electron- 
phonon (e-ph) interactions in metals and doped insula- 
tors if an adequate theory is in place. 

At present detailed experimental data on relaxation 
processes is collected for metals G, 1, B B B and high- 
temperature superconductors d, hTH) M, llO, [HI) ■ 
The pump-probe experiments are routinely analyzed in 
the framework of the so-called two temperature model 
(TTM) [H, [ii]. The model is based on the assump- 
tion that electrons and phonons are in a thermal quasi- 
equilibrium (QE) with two different time-dependent tem- 
peratures Te{t) and T/(t), respectively. The comprehen- 
sive analysis of experimental data collected at room tem- 
perature [H allowed for a determination of the electron- 
phonon (e-ph) coupling constant A of many metals and 
low-temperature superconductors in the framework of 
TTM. 

Similar experiments and their analysis were performed 
on high-temperature superconductors. The femtosecond 
time-resolved measurements on the high-T^ supercon- 
ductors TlzBazCasCuaOio [^ and YBazCusOT-s [3, [3 
found a relaxation process below T^, which is distinct 
from the equilibration of hot carriers in the normal state. 
A relatively strong e-ph coupling, A ~ 0.9 [1], a rapid 
decrease of the photo-response decay rate with decreas- 
ing temperature [l^ were found in YBa2Cu306.5, and 
the phonon bottleneck [13, [IE [13, [111 ^ biparticle 
recombination [I^, [I^ were observed below Tc- More 
recently a time-resolved photoemission spectroscopy 2C| 
and the standard pump-probe optical measurements [21] 
have been performed on Bi2Sr2CaCu208+5. Their TTM 
analysis has led to a rather weak e-ph couphng, A < 0.25. 

The pump-probe techniques have a potential to resolve 



a controversial issue on weather the e^h interaction is 
crucial [22] or weak and inessential [2^ for the mech- 
anism of high-temperature superconductivity. The pio- 
neering work by Kaganov, Lifshitz and Tanatarov [l3| 
and subsequent TTM studies are based on the assump- 
tion that electrons are in the thermal QE-state because 
the e-e relaxation time is supposed to be much shorter 
than the e-ph relaxation time. This assumption is of 
course incorrect on a femtosecond scale comparable with 
the e-e scattering time of highly excited electrons, but 
the expectation has been that deviation from QE may 
not in fact have much influence on the electron energy 
relaxation on a larger time scale [3] (for discussions of 
TTM with respect to some experiments see for example 
[11,^1). 

However later on it has been realized that nonthermal 
effects are essential even on a picosecond scale, compara- 
ble with the e-ph relaxation time, when conditions of low 
laser excitation power and relatively low temperature are 
chosen [2^. Under these conditions, the e-e collision rate 
becomes strongly suppressed as a result of the Pauli ex- 
clusion principle. Numerically integrating the Boltzmann 
equation with e-e and e-ph collision integrals Groeneveld, 
Sprik, and Lagendijk [26| have shown that the electron 
gas cannot attain a thermal distribution by e-e collisions 
on the time scale of the e-ph energy relaxation. A de- 
parture from QE leads to an increase of the e-ph energy 
relaxation time with respect to the QE expectation. As 
a consequence of this departure one might underestimate 
the e-ph coupling using TTM. 

While numerical integrations of the Boltzmann equa- 
tion can describe the time evolution of the electron dis- 
tribution function on any time scale, they require a num- 
ber of input parameters, which might be unknown a pri- 
ori. Here an analytical approach to this long-standing 
problem is developed. We reduce the integral Boltzmann 
equation to a differential Schrodinger-like equation us- 
ing an auxiliary space reciprocal to energy and find ex- 
act analytical time-dependent distributions of electrons 
in both cases of electron-electron and electron-phonon 
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relaxation. We also derive long-time relaxation rates of 
response functions and compare them with TTM. 



Here /k = df\^{t)/dt, V"c(q) is the matrix element of the 
electron-electron scattering (pseudo)potential, is the 
electron energy with respect to the equilibrium chemical 
potential. For transparency we drop the time argument 
in the distribution function. If the distribution function 
depends only on energy and time, /k = one can aver- 
age this equation over the angles of k as 



(2) 



where N{S) is the density of states (DOS) per spin, with 
the following result 

h = j j j d^'dede'K{^,^',e,e')5{i + e~^' ~e')x 

[/«'/e'(i - /«)(! - /.) - km - - , (3) 



where 



27r 



k,p,q 



We restrict our theory to relaxations involving non- 
equilibrium electron-hole excitations with energies much 
less than the equilibrium Fermi energy, Ep. Since 
the kernel fsr(^, e, e') has variation on a scale of the 
Fermi energy, one can approximate it by a constant, 
e, e') ~ K. This constant is related to the 
Coulomb pseudo-potential /ic = VciV(O), important in 
the theory of superconductivity, K « Trfi'^/2hEp. As- 
suming a low laser excitation power we linearize Eq. Q by 
introducing a small non-equilibrium correction, (/)(^, t) <^ 
1, to the equilibrium distribution, n^. 



(5) 



where — (e^' ^ + 1) ■ Keeping terms linear in 
(/)(^, t) and measuring energies in units of ksT, which is 
the only relevant energy scale of the problem, one obtains 



(6) 



II. ELECTRON-ELECTRON RELAXATION 

Let us first consider a nonthermal relaxation of the 
electron distribution function fk{t) caused by electron- 
electron collisions, which is described by the following 
Boltzmann equation, 

X (1 - 7k)(l - /p) - /k/p(l - W(l - /p-q)]-(l) 



Performing simple integrations in linearized Eq.([6l) 
yields 



cosh(^72) 



re(e) cosh(e/2) 



d^'4>{e,t)x 



sinh(£^) 2sinh(i±S^) 



(7) 



(8) 



is the familiar lifetime of electron-hole excitations in 
the Fermi liquid. Here we have used the integral 

dz\n{z)/[{z^a){z+b)] = [n^-\n{a)+ln'{b)]/2{a+b) 
with a, 6 > 0. 

The second term on the right-hand side of Eq.([7]) 
describes a source of quasi-particles due to inelastic 
electron-electron collisions. Collisions in cold degenerate 
plasmas differ essentially from quasi-elastic collisions in 
classical (hot) plasmas. In the latter the energy transfer 
is small compared with the electron energy due to a long- 
range character of the Coulomb potential, so that one can 
approximate the Boltzmann collision integral by the dif- 
ferential Landau-Fokker-Plank (LFP) equation (see, for 
example Ref . p7| ) . As one can see from Eq.Q it is not 
the case in metals. The collision energy transfer in metals 
is about the same as the excitation energy itself, which 
makes the differential LFP approximation unacceptable 
here. 

Remarkably the electron-electron collision integral ac- 
quires a differential form in a reciprocal auxiliary-time 
space introduced via the Fourier transform of Eq.(I7]), 
rather than in the energy space as in the LFP case. Let us 
consider non-equilibrium states conserving the electron- 
hole symmetry, so that the non-equilibrium part of the 
distribution is an odd function of energy, ((){—£,, t) ~ 
—<j){S,, t). If one determines a function. 



X(e,i) =0(e,i) cosh(e/2), 
then the Boltzmann equation is simplified as 



(9) 



K{kBTY 



dCxi^'^t) 



sinh(£^)' 
(10) 

We shall see below that a "bound state" of the effec- 
tive " Schrodinger" equation for the Fourier transform of 
xi^jt) corresponds to the stationary quasi-equilibrium 
distribution. Taking the Fourier transform of Eq. (|10p , we 
arrive at an exact differential counterpart of the Boltz- 
mann equation. 
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cosh (a;) 



1 



i:{x,t), (11) 
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where 



and 



(12) 



(13) 



Here another integral J^^dzz exp{ixz)/ smh{z/2) = 

27r^/ cosh^(7rx) has been used. 

The solution of Eq. (fTT|) is found as a superposition of 
normalized eigenstates, ipkix), of a textbook Hamiltonian 



k 



(14) 



where coefficients Ck are determined by the initial non- 
equilibrium distribution function (/)(^, 0) at t = 0, 



Cfe = 



dx 



de^^(x)cosh(^/2)0(^,O)e*^«/^ (15) 



The eigenstates, ipkix), and the eigenvalues, E = — fc^, 
are found from the Schrodinger equation 



dx"^ 



mirn + 1) 



cosh {x) 



V'fc(x) = fc2^fc(x), (16) 



with m = 2. This equation has a finite number of discrete 
bound states with real fc's psl ) and continuum extended 
states with imaginary fc's, k = ip {p is real ). If m is an 
integer there are m bound states with A; = 1, 2, m and 
both bound and unbound eigenstates can be expressed 
in terms of elementary functions (29j , 



ll)k{x) = AkD,nD„ 



kx 

.Die , 



(17) 



where is the normalizing amplitude and Dm — d/dx — 
TOtanh(a;). In our case (m = 2) there are two bound 
states, the even ground state with k = 2 [E = —4) and 
the odd excited state with k = \{E = —1). For relax- 
ations conserving the electron-hole symmetry the ground 
state contribution to the superposition, Ea. p4p . is inte- 
grated to zero because the initial non-equilibrium distri- 
bution is odd. On the contrary, the excited odd state 
with '4>i{x) oc sinh(a;)/ cosh^(a;) is the only state, which 
survives in Eg. p^ at i ^ cxd, so that 



V'(a;, oo) oc 



sinh(x) 



cosh^(a;) 
and (using Eq.llSl) and Eq.©) 



(a;, oo) oc 



cosh'(e/2)' 



(18) 



(19) 



which is precisely the result of the QE approximation. 
Indeed expanding the QE distribution function, /qe = 
[exp{{E - £:f )/?;) + l]"Mn powers of - T, we obtain 



4r J cosh'(C/2) 



(20) 
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FIG. 1: Relaxation of the total number of non-equilibrium 
excitations n normalized by the quasi-equilibrium number, 
uqe, for different width F of the initial non-equilibrium dis- 
tribution. 



with the same non-equilibrium correction, uqe 



^/cosh^(C/2), as in Eq.H 

The exact solution, Eq. p^ . allows us to trace 
the relaxation at any energy and at any time 
scale. In particular Fig.l represents the time evo- 
lution of the total number of non-equilibrium exci- 
tations (electrons plus holes), n = 2 Jp°° c?^0(^, t), 
for initial non-equilibrium distributions of the shape, 
0(^,0) = (7ri/2£;^/2F3)^exp[-(^/F)2], where the distri- 
bution width r is varied, but the total energy, Eq, is 
unchanged. For wide sources with large F > 1 the num- 
ber of excitations increases with time conserving the total 
energy in the process of their cooling. On the other hand, 
when most excitations at < = are created with the en- 
ergy less than fc^T (i.e. F < 1), their number decreases 
with time since their individual energies increase due to 
collisions with equilibrium electrons. One can trace the 
relaxation for any initial distribution including the case 
when a single photon initially creates a single electron 
hole pair, which cascades eventually into a large number 
of pairs until they get lost in the background thermal 
distribution. 

The asymptotic behavior of response functions can be 
readily obtained from Eq. (fT4)) by taking into account 
that only the excited bound state and the extended 
states with small p ~ {t/Te)~^/^ <^ 1 contribute to 
the sum in Eq.lTlD, when t/T^ > {l,l/T^). Substitut- 
ing the extended eigenfunctions, ippix) oc [3tanh^(a;) — 
Sip tanh(x) —p^ — 1] exp(ipx) into Eq. with Cp oc p at 
small p and integrating over p yield 



Tp{x,t) -'0(a;,oo) oc exp (-^ 



At 



(21) 



in the saddle-point approximation. Performing the 
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Fourier transform of Eq. ((2T|) with respect to x we find 



(22) 



The same result is obtained by using a r-approximation 
for the Boltzmann equation ([7|), 



0(e,t)-<^(C,oo) 



(23) 



which has the foUowing solution 

0(e, t) = <^(e, oo) + [0(^, 0) - 0(e, ^)]e-*/^=(«). (24) 



Hence one can use the r-approximation, Eq. (j24p . on the 
time scale much longer than the characteristic collision 
time, Te- However this approximation is inaccurate on 
a shorter time scale because in contrast with the exact 
solution, Eq. p^ . it does not conserve the total energy. 

Integrating Eq. yields a universal time-asymptotic 
of the total number of electron-hole excitations, 



n{t) - uqe oc 



,-t/T 
t 



(25) 



as also seen from Fig.l. 

Importantly the characteristic e-e relaxation time is 
quite long due to the Pauli exclusion principle. Using 
realistic /ic = 1 and Ep = 10 cV we estimate Tg w 1 .2 p s 
at the room temperature T = 300 K (see also Ref . |26| ) . 
which increases further as 1/T^ with cooling. 



III. ELECTRON-PHONON RELAXATION 

Now let us consider the electron-phonon relaxation de- 
scribed by the c-ph collision integral, 



2tt 



{[/k-q(l - h)N^ - /k(l - /k-q)(A^q + 1)] 
X (^(Ck - ^k-q - hu!q) + 

[/k+q(l - /k)(iVq + 1) - Ml - /k+q)7Vq] 
X (5(ek-ek+q + ^q)}, (26) 

where M(q) is the matrix element of the deformation 
potential and Nq is the distribution function of phonons 
with the frequency Wq. 

As in the former case of the e-e collisions we average 
this equation over the momentum angles using Eq.([21) 
and conventional units: 



+ - e' + - k)Nu. + h'{l - /c)]}(27) 



Here 



= ^r^5]M2(q)5(a-q-e')'5(a-e)'5(c.q-a 

(28) 



is the e-ph spectral function [l^, which has lo variation 
on a scale of the maximum phonon frequency ujo but ^ 
and ^' variation only on a much larger energy scale of the 
order oi Ep. 

Linearizing Eq. ([27|) with the help of Eq. (O yields 



0(e,t) = 27r J dojj dCQ{oJ,^,e)x 

me,t)iN^ + n^) - 0(e, t)iN^ + n_e)]5(C - ^' - hiu) + 

(29) 

Characteristic electron energies in Eq. P5)) are much less 
than the Fermi energy, so that 



(30) 



is the familiar Eliashberg function. We also assume 
that phonons are in the thermal equilibrium, N^^ = 
[exp^huj/kBT) — 1]~^, due to their fast thermalization 
caused by anharmonic interactions (i.e. phonon-phonon 
collisions) and/or due to a small size of the sample and 
the pump-laser spot allowing for a fast escape of non- 
equilibrium phonons. If this condition is not satisfied, one 
has to solve an equation for the non-equilibrium phonon 
distribution coupled with Eq. iP^]) . which is outside the 
scope of this paper. Under these assumptions Eq.f^H) is 
transformed into a form similar to the e-e collision inte- 
gral Eq.®, 



where 
1 



TphiO ftcosh(C/2) 

cosh(^'/2) 



d^stgni^ - e) 



2 sinh( 



^^0(^',t),(31) 



2-KkBT 



dwa^F 



ksTuj 



sinh^(|)tanh(f ; 



sinh(f)cosh(f) cosh(^i^)cosh(^^) 

(32) 

and energies are now measured in units of ksT. The 
Eliashberg function is quite complicated in real metallic 
compounds because of their complex lattice structures. 
This complexity can be avoided in a high-temperature 
regime, k^T ^ fkuu and in an opposite low-temperature 
regime, fcsT <C huo- 



A. High-temperature electron-phonon relaixation 

As shown by Allen [l3| the energy relaxation in TTM 
has a particularly simple form in terms of the moments 
of a'^F{uj), 



A(u;") = 2 duo- 
Jo 



'-F{uj)u 



(33) 
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where the coupHng constant A, which determines the crit- 
ical temperature in the BCS superconductors, is 



A 



dhJ . 



(34) 



Here we derive a high-temperature LFP-type equation 
for the non-equihbrium part of the distribution function 
0(^, t) to compare our exact approach with the TTM 
results [mil. 

At high temperatures the Eliashberg function is a nar- 
row function on the temperature scale, so that one can 
apply a quasi-elastic approximation expanding the e-ph 
collision integral, Eq. (P5|) or Eg. !}?!]) , in powers of the 
phonon energy, hw <C ksT . The zero-order elastic terms 
are canceled out because the distribution function de- 
pends on energy only, while the next order terms yield 
the LFP-type differential equation. 



7 



d_ 



tanh(e/2)0(O 



d_ 



(35) 



where 7 — 7r?iA(w^)/fcsT. Apart from a numerical coef- 
ficient of the order of 1 the characteristic e-ph relaxation 
rate 7/2 is about the same as the TTM energy relaxation 
rate jt = ShX{u!'^)/'KkBT, at high temperatures. 
Indeed multiplying Eg. ([55)) by and integrating over all 
energies yield the rate of excitation energy relaxation. 



Ee{t)^~l detanh(e/2)0(e,t). 



(36) 



where Ee{t) — X^<iff0(O- If replace tanh(^/2) in 
this equation by its argument assuming that 0(^, t) has 
its characteristic energy width of the order of 1, then 
Ee{t) « —{'-f/2)Ee- Hence the excitation energy relaxes 
as Ee{t) oc exp(— 7^/2) almost independent on a partic- 
ular shape of the non-equilibrium distribution. 

To verify the numerical coefficient, one can substi- 
tute the TTM distribution (pTTMi^t) = [{Te{t)) - 
r)/4T]^/cosh2(C/2) (Eq.dini)) into Eq.® to convert 
this into the temperature relaxation rate: 



Te{t) 



1 



l{Te - T) 



dxx tanh(x) / cosh^ (x) 
dxx"^ / cosh^(a;) 



(37) 



Eq. ljST]) is precisely the same as the TTM temperature 
rate [11], Te{t) = -^T{Te - T), since the ratio of two 
integrals in Eq.dST]) is G/tt^. 

According to Eq. (|55|) deviation of 0(^,t) from quasi- 
equilibrium population does not have much influence on 
the energy relaxation. Hence TTM [l^ is the ade- 
quate approximation at high temperatures, which agrees 
well with experimental observations in conventional met- 
als where Debye temperatures are rather low 



B. Low-temperature electron-phonon relaxation in 
poor metals 

Characteristic phonon frequencies are exceptionally 
high in many advanced materials like copper oxides. 



^ri/fcs ^ 400 1000 K, so that the low-temperature 
regime, fc^T ^ hwo is of great importance. Since 
all dimensionless energies in Eg. ([51]) are of the order 
of unity one can apply a low-frequency asymptotic of 
a'^F{uj) — \n[uj I lodY^ 12 in this regime. The exponent n 
depends on impurities, disorder, and sample dimensions: 
n = 2 in clean bulk crystals while n = 1 in disordered 
metals due to a phonon damping [sol . [sij and in metal- 
lic films [HI. Then Eg. l|3T|) can be Fourier-transformed 
into the Schrodinger equation using the Fourier trans- 
form Eq.lini) of x(C, t) = t) cosh(f /2). 

In the poor-metal case n = 1 the equation for x(^, t) is 
almost the same as in the e-e case, Eg.dTU]). apart from a 
numerical coefficient in front of the integral term. 



TTX{kBT)^/h^UJD 

1 



sinh(S^)' 



(38) 



where the integrals dww/ sinh(ti;/2) cosh(aj/2) = 
7r^/2 and Jg°° (iijja;tanh([jj/2)/cosh[(Li; -I- (^)/2] cosh[(a; — 
^)/2] = ^2/2sinh^(^/2) have been used. The difference 
in the numerical coefficients in front of the integral terms 
originates in different statistics of scatterers, which are 
bosons in the e-ph case and fermions in the e-e case. At 
low temperatures the e-ph relaxation time has the same 
energy and temperature dependence as the e-e relaxation 
time Eq.dH), 



(39) 



We also notice that the temperature dependence of the e- 
ph relaxation rate at low temperatures, l/Te-ph oc is 
qualitatively different from its temperature dependence 
at high temperatures 7 oc l/T. 

The Fourier transform of Eq. (l38|) yields the 
Schrodingcr-likc equation 



Tph1p{x,t) 

where 



92 



dx"^ cosh^(a 



1 



Tph 



2h^ujD 
^^XiksTY 



^{x,t), (40) 



(41) 



Different from the e-e case [Eq. p6)) with m — 2] the 
steady-state Schrodinger equation 



52 



.9x2 cosh2(a 



ipkix) = k'^^pkix) 



(42) 



has only one bound (ground) state, k ~ 1, and itinerant 
states with k — ip, 



^k{x) — Ak[k — tanh(a;)]e 



kx 



(43) 
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in the e-ph case (m = 1). Only itinerant states con- 
tribute to the superposition Eg. pi)) and determine the 
time relaxation of the distribution function because the 
contribution of the even ground state is integrated to 
zero and there is no excited odd state here. As the result 
the non-equilibrium part of the distribution function and 
the number of excitations relax with characteristic time 
Tph to zero rather than to any quasi-equilibrium state as 
shown in Fig. 3 and Fig. 4 by lower curves (a = 2). Their 
time asymptotic is found using the saddle-point approx- 
imation as in the case of the e-e collisions, 



and 



n{t) oc 



t ■ 



(44) 



(45) 



The time evolution of n(t) is widely independent of the 
width F of the initial distribution function at t = as one 
can see comparing the lowest curves in Fig. 3 and Fig. 4. 
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FIG. 2: The effective "potential" energy of the Schrodinger 
counterpart of the Boltzmann equation with the e-ph collision 
integral in clean metals compared with the potential in poor 
metals. 



Low-temperature electron-phonon relaxation in 
clean metals 



The low-frequency Eliashberg function is quadratic as 
a function of frequency, a'^F{uj) = X{(jj/lud)^ in clean 
crystalline metals, which makes an analytical expression 
for the Fourier transform of the Boltzmann equation (j3ip 
unavailable in terms of elementary functions. However 
we can approximate all relevant integrals numerically as 



duj 



sinh(a;/2) cosh((jj/2) 



..414, 



u? sinh^(^/2)tanh(tj/2) 
cosh(^)cosh(i^) 



+ 0.027r, 



and 



, cj^ cos(a;x/7r) 3tt^^^, . 
smh{uj/2) 2 ^ ' 



where V(x) is shown in Fig. 2. Then the corresponding 
Schrodinger-type equation for the Fourier transform of 
x{£,, t) becomes 



where now 



^ + 0.178^ + 0.568 



'ph 



37r3A(fcBr)3' 



(46) 
(47) 



The effective "potential" energy differs only marginally 
from the poor metal case. Fig. 2. At large t corresponding 



to large x in Eq. the forth derivative of the low-energy 
extended eigenstates is small. Hence the asymptotic be- 
havior of response functions in clean metals is qualita- 
tively about the same as in poor metals. 



n{t) oc 



exp(-0.568i/r,l) 



t 



(48) 



but the temperature dependence of the e-ph relaxation 
time is more pronounced, r^^ oc l/T^. in principle, the 
clean-metal "potential". Fig. 2, could have "resonances", 
states that are in the continuum but take a long time to 
leak out resulting in some quantitative differences with 
the poor-metal relaxation. 



IV. LOW-TEMPERATURE 
ELECTRON-PHONON RELAXATION 
COMBINED WITH ELECTRON-ELECTRON 
RELAXATION 

Finally let us combine both collision integrals into one 
Boltzmann equation. Performing its Fourier transforma- 
tion as described above in Sections (|TT1 yields the fol- 
lowing Schrodinger-like equation in the poor-metal case: 



Tli'ix,t) 



where 



and 
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dx^ cosh^(a;) 



1 



ip{x,t), 



a 



Te + Tph 



2Te -I- QTph 
Te + Tph 



(49) 



(50) 



(51) 
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FIG. 3: Relaxation of the total number of non-equilibrium 
excitations n for different e-e and e-ph scattering times char- 
acterized by parameter a, Eq. (|5H) . for the width F = 10 of 
the initial non-equilibrium distribution. 



There are two bound states of the corresponding 
steady-state Schrodinger-like equation 



52 



dx^ cosh2(x) 



-0(2;) = —Eip{x) 



(52) 



because a is larger than 2 but smaller than 6 [2g]. The 
excited odd state has the eigenfunction 



'ipiix) oc 



tanh(a;) 



[cosh(x)]l^il'^' 



and the energy 



E, = — 



(l + 4a)i/2 



(53) 



(54) 



which determines the asymptotic behavior of all linear re- 
sponse functions. In particular the number of excitations 
decays at large t as 



n{t) cx exp 



il + Ei)t 



(55) 



When both relaxations are involved the time evolution 
of n(t) calculated using Eq.([49|) with different initial dis- 
tributions differs qualitatively from TTM relaxation as 
shown in Fig. 3 and Fig. 4. In fact electrons cannot attain 
the thermal quasi-equilibrium at any a less than 6 in 
agreement with the numerical results of Ref . [5^ . More- 
over the exact relaxation rate j = 1 + Ei depends on the 
ratio of the electron-electron relaxation time, Eg. p^ . 
and the electron-phonon relaxation time, Ea. (|41|) . 



Te _ 2A Ep 

Tph 71"^? fi^n ' 



(56) 



Using Eq. (|55|) we find 



7 = c(r) 



where 



c(r) 



3^/(1 + r)(25 + 9r) - 7r - 15 



2r 



(57) 



(58) 



This coefficient changes from c(r) = 1 at r = 00 up to 
c(r) = 8/5 at r = 0. 

The TTM relaxation rate jtiow at low temperatures is 
readily obtained with the Eliashberg function a'^F{(jj) = 
XiLj/2LJD using Eqs . (4,10) of Ref. p^. \3^ . Linearizing 



Eq.(lO) of Ref. [lj| with respect to the temperature dif- 



ference Te{t)-T <^T yields 

47r3A(fcsT)2 



iTlo 



UJd 



(59) 



The ratio of our exact relaxation rate to the TTM rate 
is 



7 



iTlow 



5 ( ^ 



(60) 



If e-e collisions are much faster than e-ph collisions (r — 
0), this ratio is 1, justifying the TTM approximation also 
at low temperatures in the limit t — > cx), r — *■ . However 
at low temperatures r is not necessarily small as assumed 
in TTM even at small \ <^ jic because the Fermi energy 
in Eq. ()56|) is often much larger than the phonon energy. 
Just the opposite limit ?' — > 00 is feasible at a sizable A. In 
this limit the exact relaxation rate is slower than the low- 
temperature TTM rate, ^ /^tiow — 5/8, so that one may 
underestimate the electron-phonon coupling constant by 
about two times using TTM. Also an illegitimate fitting 
of experimental rates measured at temperatures below 
hwD/ks with the theoretical high-temperature TTM rate 
7t uM (s^^ section HII AI) may underestimate A by about 
{hwu/T^ksT)'^ times in poor metals and much more in 
clean metals. 



V. CONCLUSIONS 

In conclusion, using the auxiliary Fourier transform we 
have mapped the linearized Botzmann equation with the 
electron-electron collision integral onto a Schrodinger-like 
equation with imaginary time allowing for a simple ana- 
lytical solution. A similar mapping is also found for the 
electron-phonon collision integral at low-temperatures 
both in poor and clean metals. We have analyti- 
cally traced the time and energy evolution of the non- 
equilibrium electron distribution function on any time 
scale and found its asymptotic relaxation rate at t — *■ 00. 

A low-temperature relaxation rate strongly depends 
on the temperature: 7 oc and 7 oc in poor and 
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FIG. 4: Relaxation of the total number of non-equilibrium 
excitations n for different e-e and e-ph scattering times char- 
acterized by parameter a, Ea. (|51|l . for the width F = 2/3 of 
the initial non-equilibrium distribution. 



clean metals, respectively. The Pauli exclusion princi- 
ple slows down e-e relaxation, so that e-e and e-ph col- 
lisions are strongly entangled at low temperatures. We 
have shown that electron gas cannot attain a thermal 
quasi-equilibrium distribution by e-e colhsions, Figs. 3, 4, 
in agreement with earlier numerical integrations of the 
Boltzmann equation (26j . The rate of return to the equi- 
librium is not governed solely by electron-phonon pro- 
cesses, but also involves the electron-electron relaxation 
time, 7 — c{r)jTiow, via the coefficient c(r) which de- 
pends on the ratio r of the e-e collision time to the e-ph 
collision time. The exact relaxation rate 7 recovers its 



quasi-equilibrium TTM value 7 = JtIow only in the limit 
of the negligible e-ph coupling, A <^ it nlhw d / '2-E f ^ 0.01 
. In poor metals the physically realistic ratio r is large at 
low temperatures and the exact relaxation rate is slower 
than the TTM rate, 7 = 57t/o«,/8. 

At high temperatures, T 3> hLUu/ks, we have reduced 
the e-ph collision integral to the differential Landau- 
Fokker-Plank form. Using this form we have shown that 
the deviation of the electron distribution from quasi- 
equilibrium population does not have much influence on 
the energy relaxation, so that TTM [H [14,] is a reliable 
approximation at high temperatures. 

Our theory opens up a perspective of determinations 
of both important microscopic parameters A and /Xc us- 
ing single-parameter (a) fitting of response functions in 
pump-probe experiments at low temperatures. Figs. 3,4. 
It also allows for an analytical approach to the integral 
Boltzmann equation for the case of a steady-state source 
of excitations as in a current-carrying state. In the lat- 
ter case relaxation times could be different because the 
current carrying state does not have a distribution func- 
tion that depends only on energy, as assumed here. The 
theory could be further extended beyond the assump- 
tion that phonons remain in equilibrium by including a 
linearised Boltzmann equation for the non-equilibrium 
phonon distribution function. 
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